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We study the gravitational collapse of an inhomogeneous scalar field with quantum gravity cor- 
rections associated with singularity avoidance. Numerical simulations indicate that there is critical 
behaviour at the onset of black hole formation as in the classical theory, but with the difference that 
black holes form with a mass gap. 
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I. INTRODUCTION 

One of the outstanding problems in theoretical physics 
is the incomplete understanding at the quantum level of 
the formation, and subsequent evolution of black holes 
in a quantum theory of gravity. Although a subject of 
study for over three decades, it is fair to say that, in 
spite of partial results in string theory and loop quan- 
tum gravity, there is no widely accepted answer to many 
of the puzzles of black hole physics. This is largely be- 
cause there has been no study of quantum dynamical col- 
lapse in these approaches. Rather, progress has focused 
mainly on explanations of the microscopic origin of the 
entropy of static black holes from state counting P^l A 
four-dimensional spacetime picture of black hole forma- 
tion from matter collapse, and its subsequent evolution 
is not available in any approach to quantum gravity at 
the present time. 

This paper describes an attempt to address this prob- 
lem in the context of Hawking 's original derivation of 
black hole radiation: spherically symmetric gravity min- 
imally coupled to a massless scalar field. This is a non- 
linear 2d field theory describing the coupled system of 
the metric and scalar field degrees of freedom. 

Gravitational collapse in the classical the ory in this 
model has been carefully studied numerically^], but its 
full quantization has never been addressed. The classical 
results are well known; the onset of black hole formation 
is characterized by a scaling law for black hole mass and 
a self-similar behavior of the field variables. 

In semi-classical theory, Hawking 's calculation^ uses 
the eikonal approximation for the wave equation in a 
mildly dynamical background, where the dynamics cen- 
ters on the surface of a star undergoing collapse. The 
essential content of it is the extraction of the phase of 
the ingoing mode from an outgoing solution of the scalar 
wave equation as a classically collapsing star crosses 
its Schwarzschild radius. According to this calculation, 
black holes create particles that originate near the event 
horizon. The approximation breaks down at the late 
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stages of evaporation, where quantum gravity effects are 
expected to become important. 

It is likely that a complete understanding of quantum 
dynamics in the gravity-scalar field systerrP will resolve 
at least some of the outstanding problems of black hole 
physics associated with the late stages of black hole evo- 
lution. This has been the motivation for recent work on 
the development of a quantum framework for studying 
this modeP and for a proposal of a picture of the evolu- 
tion process that incorporates singularity avoidance^. 

In the following we give a summary of the results of 
classical gravitational collapse of a scalar field in spher- 
ical symmetry, together with some details of numerical 
techniques. We then make a quantum gravity motivated 
modification of the equations of the model designed to 
take into account "singularity resolution." Numerical 
simulations of these equations indicate that there is crit- 
ical behaviour at the onset of black hole formation as 
in the classical theory, but that black holes form with 
a mass gap dependent on the fundamental discreteness 
scale and initial data parameters. 

II. CLASSICAL COLLAPSE 

The Einstein equations coupled to a massless mini- 
mally coupled scalar field take the form 

Rab = 8Trd a (f)d b (j) (1) 

In spherical symmetry the resulting equations may be 
written with the metric ansatz 

ds 2 = -f(r, t) 2 dt 2 + g(r, t) 2 dr 2 + r 2 dn 2 , (2) 

with <p — <j)(r,t), or in double null coordinates as 

ds 2 = -4a(u, v)dudv + r 2 (u, v)dfl 2 . (3) 

with <p = <f>{u, v) . In either form these can be rewrit- 
ten as two first order evolution equations and one con- 
straint equation. Numerical studies ha ve u tilized both 
forms for the four-dimensional equations^!, and for the 
three-dimensional case with negative cosmological con- 
tant where collapse pro duces the Banados-Teitelboim- 
Zanelli (BTZ) black hole!±2in] 
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The main result is that gravitational collapse at the 
threshold of black hole formation exhibits critical be- 
haviour. For a parameter a in the initial data for the 
scalar field, numerical simulations give a mass formula 
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where a* is a critical value of a above which black holes 
form. Both a* and the critical exponent 7 are numerically 
determined. Furthermore the field variables exhibit a 
discrete self-similarity of the form 



0M) = 0(re A ,te A ) 
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where the constant A is also numerically computed. 

In the parametrization in eqn. ^ with a(u,v) :— 
g(u,v)r'(u,v), where / denotes the derivative with re- 
spect to v, the field equations in fo ur dimensions may 
be written in the compact for m 12 * 13 ! 
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where dot denotes partial derivative with respect to u, 
and we have defined 
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A n umerical integration scheme for these 
equation d 9 * 11 * 12 * 1 ^ proceeds by using a "space" v 
discretization 



h(u, v) — * hi(u), r(u, v) — > n(u) 
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to obtain a set of coupled ODEs, where i = 0, • • • N spec- 
ifies the v grid. Initial data for these two functions with, 
suitable boundary conditions, are prescribed on a con- 
stant u = uq initial slice, from which the functions g 
and g are constructed using a Simpson's rule integration. 
Evolution in u "time" is performed using the forth order 
Runge-Kutta method. 

For a physical picture, it is convenient to specify the 
initial scalar field configuration as 4>{r{v, ito), ito), since 
the metric function r is a dynamical variable in these 
coordinates. We use the data 



r(u ,v) = v 
(u Q ,r(u ,v)) = ar 2 exp 
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FIG. 1: A log-log plot of black hole horizon radius as a func- 
tion of initial data amplitude a for vanishing fundamental 
discreteness scale L. This is the expected classical result with 
slope 7 = 0.375. The oscillation of points abou t the best fit 
line is a known feature of classical collapseP^ 



At each u— step of the evolution using Eqns. 
function 
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is computed. Its vanishing signals the formation of an ap- 
parent horizon. For each run of the code with fixed am- 
plitude a, this function is scanned from larger to smaller 
radial values after each Runge-Kutta iteration. Evolu- 
tion is terminated if a root is detected within a threshold 
value. The corresponding radial value is taken as the ap- 
parent horizon radius. This procedure is repeated for a 
range of a values until its critical value a* is determined 
by a bisection method. For the subcritical case all ra- 
dial grid points evolve to zero without horizon detection. 
This is the signal for complete reflection of the scalar field 
pulse. 

Figure 1 displays typical results of a run of such a code 
for the Gaussian scalar field profile specified above. The 
results are in accord with other studies and serve as a 
useful check of our results. The slope of the line in this 
figure is 7 = 0.375 giving a mass formula of the type Q 
with the critical value a* = 0.0220025 (for a = 1 and 
tq = 1 in the initial data (13 1. The oscillation around 



The initial data parameters are a, a, tq. 



the best-fit line is also a known feature of the critical 
behavior associated with classical collapseP This data is 
actually produced from a code for the quantum gravity 
corrected equations (to be described below) , but with the 
fundamental discreteness scale L in these equations set 
to zero. 
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We turn now to describing an approach for obtaining 
quantum gravity modifications of the above equations 
that incorporate singularity avoidance. 



III. QUANTUM GRAVITY CORRECTIONS 

The origin of quantum corrections to gravitational col- 
lapse due to fundamental discreteness may be understood 
by considering quantum systems on a lattice. In spatial 
lattice based quantizations of field theory, local config- 
uration variables are sampled discretely and momentum 
variables are realized indirectly via translation operators. 
The lattice may be uniform (equispaced) or irregular. 
One can construct a non-separable Hilbert space whose 
basis elements are labeled by the lattice chosen, in ad- 
dition to other quantum numbers. Such a space may be 
viewed as the (infinite) sum of the Hilbert spaces, one for 
each lattice. 

The kinematical Hilbert space of loop quantum grav- 
ity (LQG) is such a space. The classical variables used in 
LQG, the holonomy of a connection and a surface observ- 
able associated with its conjugate momentum, are non- 
local. Their quantization is realized on a space with basis 
states labelled by graphs embedded in a three-manifold. 
Operators may also be labeled by graphs, and hence con- 
nect states with different graph labels. An example in 
LQG is a proposal for a "graph changing" hamiltonian 
constraint,^, which is defined to carry an edge that gets 
attached to a graph state. Related quantizations are 
available for scalar field theory,^ and a similar construc- 
tion exists for fcrmionic theories. 



A. Quantum Mechanics on a Lattice 

There are two types of quantum gravity effects that 
can arise from fundamental discreteness. These come 
from the way inverse configuration and momentum op- 
erators must be defined on a lattice. Both are readily 
demonstrated by considering the mechanics of a particle 
on a lattice. 

A feature of spatial lattice theories is that there is a 
natural way to define inverse configuration operators via 
finite differencing. For a particle moving on the real line, 
a quantization on a uniform lattice sampled from the 
real line with points x n — na (n £ • ■ ■ — 1, 0, 1 • ■ • ) has an 
associated Hilbert space \n) with inner product 



(m\n) = S mn 
on which we have the operators 



Uxip)\i 



exp(iXp)\n) = \n— 1) 



(15) 

(16) 
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We can define a densely defined operator corresponding 
to the variable or other inverse configuration op- 

erators, by realizing a lattice finite difference scheme for 




FIG. 2: The eigenvalue of the operator — l/x Eqn. (201 (dot- 
ted line), and the function f(x) = —l/x for A = 0.1. 



the derivative of f'(x n ) such as 

df(x n ) f(x n + A) - f(x n - A) 



dx 



2 A 



(18) 



on the Hilbert spaced For example applying it to the 
function f{x) = \f\x\ gives the scheme 




^ (Vk + a|- VK-ai)' 



(19) 



The right hand side of this expression is the eigenvalue 
of the operator 



1 

R 



x\u{-ul^/\x\u x 



u 
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in the basis \n). The eigenvalue is bounded, and hence 
may be viewed as realizing an aspect of "singularity 
avoidance" in a theory with a fundamental discreteness 
scale. 

The finite difference idea captures the essence of 
the definition on the inverse triad operator in lqgP 
and the inverse scale factor operator in loop quantum 
cosmology.^ In the general case, without symmetry re- 
duction, such operators are not bounded unless these is 
restriction to suitably defined semiclassical subspaces.^ 
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Furthermore, it is clear that there are as many such oper- 
ators as finite difference schemes in addition to the free- 
dom in the choice of function f(x); for instance a more 
general case is f(x) = Ix] 1 ^ for k > 1. These are all 
examples of quantization ambiguities. 

So far we have seen an example, in a simplified con- 
text, of an operator and its eigenvalue that arises from 
the necessity of defining operators corresponding to in- 
verse configuration variables. In quantum gravity there is 
a need to define operators corresponding to inverse metric 
functions such as the determinant of the 3-metric, which 
arises in the hamiltonian constraint. A second type of 
modification comes from the fact that momentum oper- 
ators do not exist on a lattice. As already noted these 
must be defined indirectly using the translation operators 
U\(p). One such definition is 

:= ^ {Ux - Ul) (2D 

Such expressions lead to "momentum corrections" due 
to fundamental discreteness, which in quantum gravity 
come from the variables canonically conjugate to the spa- 
tial metric - the ADM momentum. The momentum op- 
erators can also be used to define lattice creation and 
annihilation 

A±=x±i Px . (22) 

In the following study of gravitational collapse of a 
scalar field, we consider only the former type of correc- 
tion, with the understanding that a complete treatment 
would require both. The generalizations of the operators 
to the field theory case are similar. 



The gravitational phase space variables for the model 
are the configuration variables (R, (f>) and their conju- 
gate momenta (Pr,Pa). The basic Poisson bracket we 
consider is 

{R f ,e iXP ^} ee [J™ Rf dr,e iXP *^\ 

= i2G\f{r)e lXPR{r \ (23) 

where / is a suitable smearing function. 

This Poisson bracket may be represented on a Hilbert 
space with basis 

| e *'£* a k P R (x k ) ^ e iL 2 biP+to)^ 

= \a± .. .a Nl \bx ... b N2 ), (24) 

where the factors of L in the exponents reflect the length 
dimensions of the respective field variables, and a k , bi are 
real numbers which represent the excitations of the scalar 
quantum fields R and <fi &t the radial locations {x k } and 
{t/z}. The inner product on this basis is 

(ai . . . a Nl ; bi, . . . b N2 \a\ . . . a' Nl ; b\ . . . b' N2 ) 

— ■ ■ ■ 3b N2 ,b' N2 ; 

if the states contain the same number of sampled points, 
and is zero otherwise. 

The action of the basic operators are given by 

Rf \a% ... G&iVi) b\ . . . bN 2 ) 

= L 2 ^2a k f(x k )\ai ...a Nl ;b x . ..b N2 ), (25) 

k 



B. Gravity-Scalar Field Model 

The application of this type of procedure to the 
gravity-scalar field collapse problem in spherical sym- 
metry requires a lattice-based Hamiltonian quantization 
scheme. This has been developed recently.^ It contains a 
prescription for defining inverse field operators to repre- 
sent variables such as powers of 1/R(r,t), where R(r,t) 
is the metric variable that is the measure of the size of 
spheres, (r and t are radial and time coordinates.) Such 
factors arise in the Hamiltonian constraint and the field 
equations, as is evident from Eqn. ([6]), (where we have 
referred to this function as r(u, v)). 

The quantization route we follow is unconventional in 
that field momenta are not represented as self-adjoint op- 
erators; rather only exponentials of momenta are realized 
on the Hilbert space. This is similar to what happens in 
a lattice quantization, except that, as we see below, ev- 
ery quantum state represents a lattice sampling of field 
excitations, with all lattices allowed. Any numerical com- 
putation is of course only possible if it is restricted to a 
fixed lattice subspace. In this sense non-separable Hilbert 
spaces are not computationally relevant. 



e i*jPR(.*j)\ ai ...a Nl ;bi 
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■ ajvi ; bi 



>N 2 



(26) 



where aj is if the point Xj is not part of the original basis 
state. In this case the action creates a new excitation at 
the point Xj with value — Xj. These definitions give the 
commutator 



R fje i\P R (x) = -\f(x)L 2 e iXPli( - x h 



(27) 



Comparing this with ( 23 ) , and using the Poisson bracket 
commutator correspondence ih{ >}<-*•[)] gives L — 
y/2lp, where lp is the Planck length. There are similar 
operator definitions for the canonical pair (0, P^). 

This quantization is distinct from the LQG inspired 
"polymer" approach to field theorySS where it is the con- 
figuration variables that are represented in exponential 
form, following the representation of holonomy operators 
in LQG. The present approach is more akin to conven- 
tional quantization in spatial lattice field theory. It may 
be viewed as the "dual" of the polymer approach. 
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C. Inverse R operators 



The functionals 



R 



S 



drfR 



(28) 



used as the configuration variables serve as our starting 
point for defining inverse operators. The procedure for 
doing this is similar to that of the particle example. Clas- 
sical identities such as 



1 



R 



i\Gf(r) 



-i\P R (r) 



R 



/ 



JXP R (r) 



(29) 

where the functions / do not have zeroes, may be used 
to define operators. For example a suitable choice is 
a sharply peaked Gaussian at the point r k , which also 
serves to localize the operator. (It is useful to choose the 
Gaussians such that /(r fc ) = 1.) 

The representation for the quantum theory described 
above is such that the operator corresponding to Rf has 
a zero eigenvalue. Therefore we represent using the 

r.h.s. of (29 1. The corresponding operator 



1 



\R f \ \PpXfixj) 



e -i\P R (xj) 



R f \, e iAP «(^) 



(30) 

is densely defined and bounded. This may be illustrated 
with the basis state 



\ a o) 



Ja o P R (r=0)\ 



(31) 



which represents an excitation ao of the quantum field 
Rf at the coordinate origin: 



Rf\a ) 

T , , 



(2l 2 P )f(0)a \ao), 



(32) 



\Rf\ 



\ a o) = 



(A/p/(0)) 2 



«o| 1/2 -|«o-l| 1/2 ) loo) 



which is clearly bounded. If there is no excitation of Rf 
at the coordinate origin, i.e. ao = 0, the upper bound 
on the eigenvalue of the inverse operator is 2/\ 2 l 2 p . (In 
the units we are using, Pr is dimensionless, so A is a 
dimensionless number, which can be taken to be unity.) 

A symmetrical version of this operator is defined as for 
the particle case in Eqn. (20 1. Its eigenvalue on the basis 



state above is of the same functional form as in Eqn. ( 19 ) 
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D. Semiclassical states 

The basis states and their associated eigenvalues are 
not the ones we consider useful for quantum corrections 



to the collapse problem discussed in Sec. 2. Rather we 
would like to find suitable semiclassical states 



\R(r k )) 6 



(34) 



where R(r k ) is the sampling on a uniform radial lattice 
of any prescribed function R(r) , with the properties that 
for each lattice point rjt we have 



Rf) =R(rk) 



(35) 



and 



\Rf\ 



(A?P./(r fe )) 2 

x (y/\R(r k ) + Aj - y/\R(r k ) - A|)(36) 



The right hand side of the latter function has the prop- 
erty that for R(rk) >> A, it behaves like 1/R(r k ) but 
has a different functional form for R(r k ) ~ A, just as in 
the particle mechanics case. 

It is possible to explicitly construct such statesPD The 
idea is to associate semiclassical states for a single par- 
ticle on a lineal with the lattice points {rfc}, such that 
a state at point n is peaked at coordinate value R{fi). 
The state \R(rk)) is then defined to be the product of 
such single particle states. It represents a quantum state 
corresponding to the classical profile R(r). (It is conve- 
nient to take f(rk) — 1 and set Xl p — L, which we do in 
the following. ) 



E. Modified Collapse Equations 

A derivation of quantum gravity corrections to the 
field equations requires a number of inputs depending 
on the approach taken. One approach is to make use 
of semiclassical states peaked on classical configurations. 
Given a definition of constraint operators on a kinemati- 
cal Hilbert space, one then computes expectation values 
of the constraints in such states. This results in quan- 
tum corrected "effective constraints," which to leading 
order in H are the classical constraints. The idea is to use 
these new constraints to derive evolution equations. For 
example, if H(x,p) is such an operator, its expectation 
value in a semiclassical state \x,p) would give a func- 
tion H(x,p). This function would then be used to derive 
quantum corrected equations. This is the approach we 
take. Similar methods has been used to derive effective 
equations in loop quantum cosmology,^ and partly form 
the motivation for our work. 

It is apparent from sample computations that this ap- 
proach gives the classical equations in suitable limits. For 
example if there is a 1/x factor in a function H(x,p) (as 
is the case in the Hamiltonian constraint), the equation 
of motion also contains a power of 1/x. The correspond- 
ing quantum corrected equation will have a power of the 



6 



spectrum function of the corresponding operator as in 
Eqn. (19 1. Given the form of this function, it is appar- 



ent that the large x limit gives the classical equation as 
shown in Figure 2, and that the small x region represents 
a "repulsion," if the eigenvalue represents a potential en- 
ergy. 

A possible (and perhaps obvious) criticism of this ap- 
proach to obtaining quantum corrected equations is that 
one is deviating from the ideal of solving the quantum 
constraints a'la Dirac and obtaining the physical Hilbert 
space before proceeding to study physical questions. This 
approach has however not yielded any results beyond 
simple mini-superspace reductions, in the metric or the 
connection variable constructions. In order to proceed it 
is important to see what quantum gravity effects can be 
extracted from a kinematical construction of constraint 
operators. The approach outlined above is one example 
of such a procedure, but it is an approximation which 
requires careful scrutiny especially with regard to the 
consistency of the quantum corrected constraints to the 
appropriate orderPA 

The eigenvalue of the 1 /R field operator on basis states 
contains the modulus of the eigenvalue of the radial field 
operator R. This causes a numerical problem in the 
quantum modified equations where it is necessary to com- 
pute derivatives with respect to a radial coordinate. To 
avoid this issue we replace the expectation value ( 36 1 by 
the smoothed version 




-(R/Lf 



(37) 



for numerical calculations. This form has the same qual- 
itative features as the actual eigenvalue, i.e. the large 
R behaviour is 1/R with repulsion at small R, as in 
Fig. 2. It is possible that an exact computation of this 
expectation value in semiclassical states will provide a 
smooth function, so that such an approximation would 
become unnecessary. However we expect that the quali- 
tative features of the result described below will remain 
unchanged. 

Modified collapse equations may be obtained by first 
computing effective constraints and obtaining equations 
of motion from them, or by directly replacing the ex- 
pectation value ( 37 1 in place of factors of 1 jr into the 
equations from Sec. II. These two procedures of course 
give different effective equations. The procedure we used 
amounts to the latter since all the corrections from the 
former approach have not been incorporated into the 
code used. The modification is such that for large spheres 
(i.e. large R), the equations converge to the classical 
equations, with quantum gravity corrections confined to 
smaller R values determined by the fundamental discrete- 
ness scale L. As already noted, we do not include the 
momentum corrections in the present work; this would 
amount to replacing momenta in the constraints such as 
Pr by the compactified form sin(APij)/Ap^ 



IV. RESULTS 

The numerical procedure used for the quantum gravity 
corrected equations is identical to the one outlined in 
Sec. II. The simulations were performed for a range of L 
values, with the integration lattice chosen such that time 
and space steps ranged from 10 -3 to 10 -4 . The initial 
data used was the same as that for the classical case. 

The results for horizon radius Rbh as a function of the 
initial data parameter a are given in Fig. 3, for funda- 
mental discreteness scale L values ranging from to 0.2. 
There are a number of points to note here: (i) L = 
gives no mass gap, and the expected classical result, (ii) 
Mass gaps at the onset of black hole formation are evi- 
dent, with the critical value of the amplitude a* depen- 
dent on L; the gaps increase with L. (iv) For amplitudes 
a >> a* the points corresponding to different L values 
begin to merge, i.e. for sufficiently massive initial data, 
the black hole masses become independent of L. This is 
expected intuitively because bigger masses should be less 
affected by smaller scale quantum effects than smaller 
masses, (v) The critical amplitude a* depends on L in 
a rather unusual way: for non-zero L, a* decreases with 
increasing L. This means that black holes form more 
readily for larger L but the effect is rather small as far 
as the parameter ranges are concerned. 

Our results in Fig. 3 can be summarized in the black 
hole mass formula 



AI 



BH 



mo(L, a) + k[a — a*(L)] 



7 (L,a) 



(38) 



in the supercritical region a > a* , where too is the mass 
gap and k and 7 are numerically determined constants. It 
is apparent from the data that the mass gap to increases 
with L, and that the exponent 7 is also a mild function 
of L. Dependance on the amplitude a is also apparent in 
that for large values, the points appear to converge. 



V. SUMMARY AND DISCUSSION 

We have described a first exploration of possible quan- 
tum gravity corrections in the gravitational collapse of an 
inhomogeneous scalar field in spherical symmetry. Nu- 
merical simulations indicate that there is critical behav- 
ior at the onset of black hole formation, and that black 
holes form with a mass gap unlike in classical theory. 

There are general grounds to expect the mass gap re- 
sult given that there is a fundamental discreteness scale 
in quantum gravity. This is simply that any physical ob- 
ject must have mass and size in units of this scale. It is 
gratifying that this is borne out, especially in the details 
of its dependence on both this scale and on the amplitude 
of the initial data. 

A mass gap in black hole formation has been noted 
before in a model where an exterior generalized Vaidya 
solution is patched to an interior Friedmann-Robertson- 
Walker modeP^. This is also related to singularity avoid- 
ance in that the interior incorporates the "inverse scale 
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FIG. 3: Black hole radius as a function of initial data am- 
plitude a for various values of the fundamental discreteness 
scale L. Black hole mass gap increases with L. 



factor" bound used in quantum cosmology! 16 * 17 ! We note 
that there is no critical behaviour in this model because 
there is only matter inflow, and a black hole always forms 
by construction regardless of the initial data. This of 
course is a feature of Oppenheimer-Snyder type models 
of gravitational collapse with pure inflow. Only models 
with both matter inflow and outflow (i.e. metric func- 
tions fully dependent on both r, t or u, v) have the pos- 
sibility of critical behaviour, with long time evolution 



leading either to black hole formation or full reflection. 

There arc a number of directions for further work based 
on the approach we have used. These include incorporat- 
ing the momentum modification into the evolution equa- 
tions, and a larger exploration of the parameter space. 
With this in place it would be especially interesting to 
explore the nature of the exactly critical solution a = a*; 
in the classical theory this is a naked singularity that re- 
alizes a finely tuned violation of the cosmic censorship 
hypothesis. With quantum gravity corrections the criti- 
cal solution might be a boson starP^l 

A further challenge is to set up a numerical procedure 
that can be used to evolve data past horizon formation. 
This will likely require coordinates such as the flat slice 
onea^l, that have the potential to reveal how the scalar 
field and horizon evolve beyond its formation. If the hori- 
zon is found to shrink, it would be accompanied by a flux 
of scalar field away from the hole region, and should leave 
a clear numerical signature in the evolution of this field. 
An application of our approach to the possibly simpler 
3— dimensional case with negative cosmological constant 
may be easier to study in this respect. 

Another direction with regard to such questions is the 
path integral. There is a derivation of black hole radiance 
in this approach^- 8 , and it would be interesting to see 
how quantum gravity corrections arise in the spin foarrPsI 
approach. Given the understanding of 3— dimensional 
gravity^ using spin foam methods, this may be a first 
example to consider. 
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